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We present a model for the present accelerating Universe and focus on the different important 
physical variables involved in the model under the phenomenological assumption A oc H 2 with a 
prescription for equation of state parameter in the form cj(t) — uio + ^Sf , where o»o and wi are two 
constants and r is a parameter having dimension of time t 2 . General expressions for the density 
parameter SI and deceleration parameter q are obtained which under specific bound reproduce some 
of the previous results. We explore physical features of these parameters which (i) provide the 
scenario of complete evolution of the cosmos with ui(t) and (ii) agree mostly with the observational 
status of the present phase of the accelerating Universe. 

PACS numbers: 04.20. Jb, 98.80.Jk, 98.80.Cq. 



I. INTRODUCTION 



Ever since the discovery of an Universal acceleration 
[E 0; g rea t efforts are going on for finding a satisfac- 
tory explanation of this effect. Most of the works re- 
lated to this problem are based on the assumption that 
an exotic type of energy, termed as dark energy, is the 
causative agent of this acceleration although alternative 
explanations are also available in the literature 
Phenomenological models are also contenders in this race 
for understanding the elusive dark energy. A number of 
phenomenological models involving time-dependent cos- 
mological term A have already been proposed and tested 
by various workers from different standpoints (for an 
overview see Overduin and Cooperstock |6j and Salmi 
and Starobinsky Q). In favor of A-decay scenario, irre- 
spective of whether they come from extended theories of 
gravity or phenomenological considerations, it is argued 
[a, Q that (i) they have been shown to address a number 
of pressing problems in cosmology; (ii) many are indepen- 
dently motivated, e.g., by dimensional arguments, or as 
limiting cases of more complicated theories; (iii) most are 
simple enough that meaningful conclusions can be drawn 
about their viability and (iv) successful implementation 
would point either towards the eventual Lagrangian for- 
mulation or a more complete theory. 



A very important parameter for the above dark en- 
ergy investigation is that of the equation of state pa- 
rameter lu. For pressureless dust, radiation, stiff- fluid 
and vacuum- fluid dominated Universe, one can assume 
the values 0, 1/3, 1 and —1 respectively for ui. Some 
other limits obtained from observational results coming 
from SN la data [|| and SN la data corroborated with 
CMBR anisotropy and galaxy clustering statistics [l(| 
are -1.67 < ui < -0.62 and -1.33 < w < -0.79 respec- 
tively. However, it is not at all obligatory to use a con- 
stant value of ui. In fact, it can be a function of time, red- 
shift or scale factor |llrll3| . But, due to lack of conclusive 
observational data which may enable us to distinguish be- 
tween a variable and constant equation of state, usually 
a constant value of u is used [H], EH ■ Coming to the case 
of particular relationship of variable ui with redshift, it 
can be linear like ui(z) — ujq+uj'z where oj' = (dui/dz) z= (] 
[_J_,[___| or non-linear as ui(z) = uiq + ui\Z / (1 + z) |18l,Tl9j. 
So, as far as the scale factor dependence of ui is con- 
cerned, the parametrization u(a) = luq + u a (l — a), where 
ujo is the present (a = 1) value and uj a is a measure of 
the time variation ui', is widely used in the literature 
[2fj| . Different forms of ui a are also available in various 
works [lj| [2l[ . Some other models such as kink model 
[22| and e-fold model 23 1 have tried to construct a mono- 
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tonic ui(a) involving four parameters which makes transi- 
tion from some asymptotic past value ui p to some asymp- 
totic future value uif. Physical supports in favor of ui a 
parametrization have been given by Linder (Tsi l24j . In 
quintessence models involving scalar fields give rise to 
time-dependent u> [25i - l2-i |. There are also some cosmo- 
logical models with a time dependent equation of state 
parameter available in the literature with Kaluza-Klein 
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metric and wormholes |3CH32| . 

Another aspect of the accelerating Universe is related 
with the signature flipping of the deceleration param- 
eter q. In fact, in the present dark energy dominated 
Universe, the second largest constituent is another dark 
component, viz. dark matter. About 25% contribution 
from the matter content of the Universe is in the form of 
dark matter. Now, in an ever accelerating Universe, hi- 
erarchical structure formation is impossible. A primary 
criteria, therefore, for the A-CDM Universe is that the 
present acceleration is a recent phenomena and was pre- 
ceded by a period of deceleration. Hence the deceleration 
parameter q must have changed its sign sometime during 
cosmic evolution. Both theoretical [33|, HH and observa- 
tional 35] works support this physical idea. So, q must 
be time-dependent so that the transition of the Universe 
from a decelerating to an accelerating phase can occur 
smoothly. This point must be kept in mind while dealing 
with dark energy investigation regarding present acceler- 
ating Universe. 

In recent years various forms of time-dependent U) have 
been used for variable A models, viz. Aa H 3 , [36l - l38| , 
A oc H [H and A oc H 2 0. Recently Ray et al. [Il[ 
have shown the equivalence of the three A models, viz. 
A oc H 2 , A oc a/a and A oc p for constant u>. It was 
mentioned in that work of Ray et al.[4l| that in a later 
work, a more accurate analysis may be made by consid- 
ering a time- varying equation of state parameter whereas 
Mukhopadhyay et al. have already investigated the 
A oc H 2 model for a particular type of variable w with the 
full physical range of a. So, the purpose of the present 
paper is to re-examine the status of the same A oc H 2 
model for the presciption oj(t) = ojq + where ujq 
and uj\ are two constants and r is a parameter having 
dimension of t 2 . The scheme of the investigation is as 
follows: In the Section 2 the Einstein field equations and 
their solutions are provided whereas Section 3 deals with 
the physical features of different parameters. In Section 
4 some comments are made on the basis of the results 
obtained in the investigations. 



II. THE EINSTEIN FIELD EQUATIONS AND 
THEIR SOLUTIONS 



The Einstein field equations are 
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where the erstwhile cosmological constant A is assumed 
as a function of time, viz. A = A(i) and the velocity of 
light c in vacuum is unity when expressed in relativistic 
units. 

In this connection it is to be mention that invariant 
property of A under Lorentz transformation is not sat- 
isfied for arbitrary systems (e.g., material systems and 
radiation). It has been argued that the energy den- 
sity of vacuum represents a scalar function of the four- 



dimensional space-time coordinates so that it satisfies the 
Lorentz symmetry [42], [4j|. On the other hand Vish- 
wakarma [44j considered a particular Ricci-symmetry un- 
der the framework of general relativity which is the con- 
tracted Ricci-collineation along the fluid flow vector and 
argues that this symmetry does demand Lambda to be a 
function of time (and space, in general). In this context 
we would also like to mention that by the application of 
renormalization group approach of quantum field theory 
in cosmology Shapiro et al. 45| have derived generalized 
energy conservation laws for variable A models. Such 
generalized laws have also been deduced by Vereschagin 
et al. based on the work of Gurzadyan and Xue [43] . 
However, comment made by Beesham [48[ and Shapiro 
et al. [45j that for simplicity, one may still retain the 
usual energy conservation law even when A is a variable. 
In the present investigation, therefore, we have used the 
same energy conservation law for solving the dynamical 
cosmological equation. 

For the spherically symmetric Friedmann-Lemaitre- 
Robertson- Walker (FLRW) metric 
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where a is the scale factor and the curvature constant k 
is — 1,0, +1 respectively for open, flat and close models 
of the Universe, the Einstein field equations ((T|) take the 
forms as 
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From equations Q and (j3| , one can arrive at the equa- 
tion El 
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where the time-dependent barotropic equation of state 
parameter u>(t) is given by uj(t) = p(t)/p(t). 

Let us use the ansatzK = 3aH 2 . This is one of the very 
well known and widely used phenomenological models 
as proposed by Carvalho et al. [Ij| and Waga[50j from 
dimensional arguments. However, Lima and Carvalho 
[5l| have considered the same model from different point 
of view. According to Vishwakarma [52j | , one can obtain 
the relationship A oc H 2 when the cosmological term A is 
expressed in terms of Planck energy density. Moreover, 
it has been demonstrated by other authors [53j, l54| that 
A oc H 2 law can also be deduced by using effective field 
theory and black hole thermodynamics. 

By the use of this ansatz A = 2>aH 2 in the equation 
(0, we arrive at the following form 
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which, on simplification becomes 



dH 
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-(l + w)dt. 



(7) 



Now, it is important to consider a suitable form of the 
equation of state parameter lo. As already mentioned in 
the Introduction that it can be time-dependent and even 
it may be a function of the red-shift z or scale factor a 



as well. For time-dependence Mukhopadhyay et al. [55 1 
suggest it in the form u)(t) = loq + u>it n which can be 
regarded as a generalization of the special form u>(t) — 
ujq + ui\t [361 ] . In the present work, let us suppose the 
structure of u as follows 



Oj(t) — 0J + 



UJ\T 



(8) 



where ujq and w\ are two constants and r is a parame- 
ter having dimension of t 2 . Here r apart from being a 
non-negative mathematical parameter, rather bears some 
deeper physical significance. It has been shown that r 
represents the time-scale of evaporation of Bose-Einstein 
condensates which include a time-dependent A [56H58l | . 
In this context it is interesting to note that Mukhopad- 
hyay et al. (5|| consider r in their equation of state in 
the form r = t/(l + w) 1 /". Thus, they argue that for 
the present dust-filled Universe (w = 0), r is equal to t 
whereas for vacuum fluid r becomes meaningless. How- 
ever, in the present case the expression for u) immediately 
suggests that for very small t, i.e. in the early Universe, 
the contributions from the second term, u>iT/t 2 , in the 
right hand side of equation © were significant. But for 
large t, i.e. for late Universe, this term tends to zero and 
consequently uj(t) converges to a nearly constant value 
ujq. This may be one of the reasons that at present we 
can not distinguish between a constant and a variable 
equation of state as mentioned in the Introduction. 

By the use of the equation © in equation ([7]) and 
solving the resulting equation we get our solution set as 
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Now, since the Universe is at present old enough, one 
can easily neglect the terms involving powers of 1/t 2 
greater than one. Thus the equations (fTU|) -(IT2l trans- 
form to 
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From equation (|T3|) it is clear that for large t, the term 
wit/(1 + oJo)t 2 is negligibly small and hence we get the 
relationship H oc Similarly from equations (|14[) and 
(|15l) neglecting the terms involving t~ 4 one can arrive 
at the relations p oc i~ 2 and A oc t~ 2 . Thus for large t 
we can recover the relationships H oc t , p oc t~ 2 and 
A oc i~ 2 obtained by Ray et al. [4l[ for constant u>. This 
proves the generality of the present work so far as the 
time variations of iJ, p and A are concerned. 



III. PHYSICAL FEATURES OF THE 
PARAMETERS 

Using equations ([T5]) - ([T5]) we can obtain the expres- 
sions for the matter-energy density f2 m and dark-energy 
density J7a as 
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Adding equations (| and ([T7]) we can obtain 



fi = SI, 
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From equations ([T5]) and ([T7]) we find that for large 
t, Sl m ~ (1 — a) and S1a — ol. For the same A model 
with constant ui, Ray et al. [4l| obtained Sl m = (1 — a) 
and SI a = ol. So, in this case also, we find that the 
present work is a generalization of the work of Ray et 
41| . Further from equation (fT5|) we get S7 ~ 1 when 
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t is very large. This result is also compatible with the 
observational results and the work of Ray et al. [4lj ]. 
Moreover, since the second term in the right hand side 
of equation (fT%j) is clearly positive, SI < 1. But, the 
sum of all types of energy densities being equal to 1, 
the term 4w 2 t 2 /(1 + wo) 2 i 4 in equation (ITS)) can be in- 
terpreted as the sum of energy densities due to radiation 
and curvature, or in absence of any curvature, the energy 
density coming from radiation only. It has already been 
mentioned that for large t, the term 4w 2 t 2 /(1 + wo) 2 £ 4 
is very small. This result is compatible with the idea 
that at present both radiation-energy density and energy 
density due to curvature are negligibly small and almost 
entire energy density of the Universe comes from matter- 
and dark-energy densities. In this context we would like 
to mention the work of Mukhopadhyay et al. [59J where, 
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FIG. 1: Variation of density parameter with cosmic time for 
the specified values of wo when Wi = 0.5 (Equation 18). 



[Q m = 0.3, ^ = 0.5] 

3 - 




FIG. 2: Variation of deceleration parameter with cosmic time 
for the specified values of luq when w m = 0.3 and uji = 0.5 
(Equation 20). 



for the A oc H 2 model with constant oj, it has been shown 
that the sum of matter- and dark energy-densities is equal 
to 1 for both early and present Universe. But this result 
Q m + fl\ = 1 is valid for the late Universe also only when 
the Universe is composed of two fluids having barotropic 
indices uj a and uib- In the present work, however, without 
resorting to any special assumption, it has been possible 
to show that with the passage of time the density pa- 
rameter fl (= Q m + smoothly moves from values less 
than 1 to the currently accepted value 1. This aspect 
is clear from the plot for f2 vs t in the Fig. 1, where 
variation of density parameter with cosmic time for the 
specified values of ojq when uj\ = 0.5 has been shown. 

Let us now calculate the deceleration parameter to see 
it's physical status. For this, by using equation (fTU|) . we 
get the expression for the deceleration parameter q as 



Q = 



3(1 - a) 



- 1. 



(19) 



From equation (H~9l we find that for small t, the term 
u>\r/t had a significant contribution and consequently 
the Universe was in a decelerating phase. But with the 
passage of time the term ojir/t 2 became negligibly small 
and hence the Universe switched over to the present ac- 
celerating epoch. This clearly shows the signature flip- 
ping of q which is essential for the A-CDM cosmology. 
Therefore, for large t one can neglect the second term 
inside the bracket in equation (|16l) . when the equation 



([T^| reduces to to following form as 
3Q 



(l +Wo ) + — 



1. 



(20) 



Now, assuming the value f2 m = 0.30. which is within 
the range f2 m = 0.330 ± 0.035 [H, [|KH52[ and cj = 0, 
the present value of the q comes out as go = —0.55. This 
value of go nicely fits with the presently accepted value of 
that parameter as go = —0.50 ± 0.05 for an accelerating 
Universe [63], |64[ (Fig. 2). This value of qo was earlier 
obtained by Ray et al [4JJ for lj = 0. Thus through the 
present work, it has been possible to recover the value of 
go as obtained by Ray et al. [4lj |. 

In the same spirit we would like to look at the equation 
of state expressed in the equation (9) and plot the vari- 
ation of equation of state parameter with cosmic time. 
The evolution of u are shown here for radiation, dust, 
and stiff fluid phase with solid, dashed and chain-dashed 
lines respectively (Fig. 3). For all these, u)\ is taken to 
be 0.5. However, it is observed that the nature of plots 
related to to does not depend on the values of u>\ . 



IV. DISCUSSION 

In the present investigation some interesting physics 
are explored by assuming a particular time-dependent 
form of the equation of state parameter u> to track down 
the time evolution of the Universe. It is shown that at 
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FIG. 3: Variation of equation of state parameter with cosmic 
time for the specified values of uio when uj\ — 0.5 (Equation 
8). 



present both radiation-energy density and energy density 
due to curvature are negligibly small such that almost the 
entire energy density of the Universe comes from matter- 
and dark-energy densities. With the passage of time the 
cosmic density parameter SI smoothly moves from values 
less than 1 to the currently accepted value 1 as shown in 
the Fig. 1. 

It has also been possible to recover the results of Ray 
ct al. 41] so far as the expressions of H, p, A, f2 m , f^A 
and q are concerned. In particular, the expression for q 
shows the required signature change during cosmic evo- 
lution proving that the present acceleration of the Uni- 
verse is a recent phenomena (Fig. 2). The calculated 
value of q is compatible with the modern accepted value 
of this parameter. The present investigation has also 



been interesting to show a glimpse of the small amount 
of radiation-energy density left over from the early Uni- 
verse in the form of CMBR and tiny fraction of energy 
density, if it exists at all at present, due to curvature of 
space. This justifies the comment made by Mukhopad- 
hyay et al. |59| that we are still not in a position to rule 
out k — ±1 cosmologies. 

It is to note that the present work has also been suc- 
cessful in justifying the idea of Ray et al [4l[ that time- 
dependent equation of state parameter is essential for 
studying the complete time-evolution of the cosmos (Fig. 
3). In connection to the work of Mukhopadhyay and Ray 
[6a ] it is to be mentioned here that for quintessence, vac- 
uum fluid and phantom energy, the rate of acceleration 
should be higher, e.g. for oj = —0.5, —1.0 and —2.0 the 
resultant deceleration parameter values are respectively, 
q = —0.775, —1.0 and —1.45. On the other hand, for 
stiff fluid (cu = 1.0), q = —0.1. This suggest that more 
smaller the value of ui the higher is the rate of acceler- 
ation resulting the so-called Big Rip 66] or Partial Rip 
[67J scenario due to divergence of scale factor. One can 
also note the work of Kuhlen et al. [68[ where they con- 
sider a range of parameter space: cj = —0.5, —0.75, —1.0, 
— 1.25 and —1.5 for their simulations to reproduce the 
time-evolution of the Universe. However, in the present 
investigation by considering an equation of state in the 
form us(t) — uiq + we are able to get the total range 
of time evolution of the Universe. Here, as can be seen 
from the Fig. 3 that starting from different values for ra- 
diation, dust and stiff fluid all the situations ultimately 
approach to the same physical feature. 

Finally, we would like to mention that in the present 
model p is positive and hence the Weak Energy Condition 
(WEC) as well as the Null Energy Condition (NEC) are 
obeyed here. 
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